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ABSTRACT 

We analyze the redshift evolution of the mass-metallicity relation in a sample of 1 10 Damped 
Lya absorbers spanning the redshift range z = 0.11 — 5.06 and find that the zero-point of 
the correlation changes significantly with redshift. The evolution is such that the zero-point 
is constant at the early phases of galaxy growth (i.e. no evolution) but then features a sharp 
break at z = 2.6 ± 0.2 with a rapid incline towards lower redshifts such that damped absorbers 
of identical masses are more metal rich at later times than earlier. The slope of this mass 
metallicity correlation evolution is 0.35 ± 0.07 dex per unit redshift. 

We compare this result to similar studies of the redshift evolution of emission selected 
galaxy samples and find a remarkable agreement with the slope of the evolution of galaxies of 
stellar mass log(M* /M Q ) « 8.5. This allows us to form an observational tie between damped 
absorbers and galaxies seen in emission. 

We use results from simulations to infer the virial mass of the dark matter halo of a typical 
DLA galaxy and find a ratio (M V i r / M*) w 30. 

We compare our results to those of several other studies that have reported strong 
transition-like events at redshifts around z = 2.5 — 2.6 and argue that all those observa- 
tions can be understood as the consequence of a transition from a situation where galaxies 
were fed more unprocessed infalling gas than they could easily consume to one where they 
suddenly become infall starved and turn to mainly processing, or re-processing, of previously 
acquired gas. 

Key words: galaxies: formation - galaxies: evolution - galaxies: high-redshift - galaxies: 
ISM - galaxies: mass-metallicity relation - quasars: absorption lines - cosmology: observa- 
tions 



1 INTRODUCTION 

The conversion of gas into stars and the ensuing enrichment of 
the primordial gas with heavier elements is a fundamental pro- 
cess of cosmic evolution, its monotonic nature making it akin to 
an arrow of time for that evolution. Heavy elements are at high 
redshifts most eas ily detected i n abso rption and Damped Lyman a 
absorbers (DLAs. lWolfe et alj dl986h ) allow us to trace the metal- 
licity of the cold gas in galaxies back to redshifts of z > 5 (e.g., 
iLu et al.ll996l ; |Prochaska et al.l2003l;lRafelski et al.l2012h . Initially 
it was thought that DLA studies directly would provide us with the 
history of metal enrichment and that models based on DLA sam- 
ples would be able to fully describe this evolution (e.g. JPei & Falj 
l l995l : ICen & Ostrikerl 19991) . However, the data displayed a surpris- 



ingly weak dependence of the DLA metallicities on redshift (for 
a review see |Pettinill200d) . Also, the population of Lyman-break 
galaxies (LBGs) were discovered and some of these were found 
to have significantly higher met allicities than foun d in the DLA 
systems at similar redshifts (e.g., IPettini et al.ll2002l) . The relation 
betwee n PL As and LBGs was subsequently the subject of several 
studies dFvnbo et al.lll999l;lBunker et alj|l999l:lwarren et al.ll200ll; 
lM0ller et al.ll2002k IWeatherlev et alj|2005l iRauch et al.ll2008l) and 
it was found that the difference between the two sets of galaxies 
can be understood from their very different selection functions (gas 
cross-section vs. UV luminosity selection). Also, it is becoming 
increasingly clear that correlations between metallicity and other 
fundamental galaxy parameters like luminosity, stellar mass, and 
star-formation rate are central to understandi ng the galaxy popula- 
tion both locally and at high redshifts (e . g..lTremonti et al]|2004l ; 
ISavaglio et all 120051 : ICalura et al.l |2009; Mann ucci et al.l 12010b . 
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There is moun ting evidence that DLA galaxies at z > 2 follow sim- 
ilar relations dM0ller et al 112004 iLedoux et al.ll2006l: iFvnbo et al 



2012) 



2008; Pontzen et al. 2008; Prochaska et all 120081 : iKrogager et al 



The Lyman break technique and searches for DLAs in quasar 
spectra represent two widely different methods for finding galaxies 
at high redshifts and as argued above they are now understood to 
select different sub-samples of the high redshift galaxy population. 
Yet, there is an overlap between the two which can be exploited to 
gain a more detailed understanding of the galaxies in the overlap- 
ping region. 

Searches for high redshift Lya emitters is a third and inde- 
pendent method which again selects an independent well defined 
sub -sample of the und erlying sample of all high redshift galax- 
ies. lNilssonetal]d2009l l201 1) found a remarkable property of Lya 
emission selected galaxies namely an apparently sudden transition 
in their dust propertie s between redshifts 2.3 and 3 (confirmed by 
ICiardulloetai1l2012l) . In a n unrelated study of highly reddened 



Fvnbo et al. ( 120131) reported a sudden drop in the detec- 



quasars 

tion rate of dusty quasar host galaxies at redshifts above 2.6 in 
good agreement with the transition of dust in Lyman a emitters 
at z ~ 2.5. 

In the same redshift range a transition in the mode of star for- 
matio n is s een to take place in gal axy formation models dOser et all 
l20ld) and lMannucci et alj (2010) reported that here a tight 3-way 
relation between mass, star formation and metallicity appears for 
the first time. This relation is thought to be caused by the inter- 
play between infall of pristine gas and the outflow of enriched gas. 
The fact that it breaks down at z > 3 suggests that the balance be- 
tween inflow and outflow may be changing above and below this 
reds hift, possibly caus ed by differences in the initial mass func- 
tion |Bailin et al. 2010 ) or the ther modynamics of the intergalactic 
med ium jD ixon & Fu rlanettoll2009h . 

ILedoux et alj d2006j) showed that there is a strong mass metal- 
licity relation for DLAs out to redshifts beyond z — 4 but they also 
reported tentative evidence that this relation was changing with red- 
shift. The aim of this paper is to further investigate and quantify if 
there is indeed a redshift evolution of the DLA mass metallicity re- 
lation, how this relates to t he redshift evolution of the mass metal- 
licity relation of galaxies dMaiolino et al.ll2008n and in particular 
how this is related to the phase transition in the physics of galaxies 
and star formation at a redshift around 2.5-2.6 as described above. 



2 REDSHIFT EVOLUTION OF THE METALLICITY 
MASS RELATION 

2.1 Sample definition 

Our basis is the sample of 70 absorbers from ILedoux et all d2006l) 
(their Table 1). In order to increase the redshift coverage over that 
sample we have searched the literature for DLAs at both very low 
and very high redshifts with published spectroscopic observations 
of sufficiently high resolution. We found a total of 14 low redshift 
DLAs (in the redshift range 0.1 - 1.5, Table 1) and we include them 
all here. At high redshifts we sele cted the 22 DLAs observed with 
HIRES from Raf elski et alj fe012h which provides a homogeneous 
high redshift sample providing good coverage back to z = 5.06 
(Table 2). In the sample we also includ ed four newly ob served 
objects from our X-shooter DLA survey dFvnbo et al since 
the parameters required were already extracted and published (Ta- 
ble 3). Our total sample therefore includes 110 systems sampling 
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Figure 1. Left frame: Metallicity versus velocity width of each system in 
our sample. The best fit lineal' relation using a slope of 1.46 (see text) is 
drawn as a red line. Right frames top to bottom: Residuals after subtraction 
of the fit with no redshift evolution (eq(2)) and after subtraction of the fit 
including redshift evolution (eq(2 and 5)). 



the redshift range z — 0.11 to 5.06. This represents an increase in 
sample size of 60% but more importantly it allows us to study the 
redshift evolution over a large fraction of the age of the universe. 

2.2 Method 

There is a relation between the velocity width of a DLA sys- 
tem and the gravita tional mass of the host dHaehnelt et alJfl998L 
IPontzen et alj 20o3) which makes it possible to use the velocity 
width, Ay, as a proxy for the mass. In this paper we shall use the 
A V as defined in Ledoux et al ]£oo3) and for consistency we care- 
fully choose suitable low ionization absorption lines as described in 
that work. To convert the published high resolution line profile data 
to electronic data we used the ADS Dexter data extraction applet 
which is a tool to extract data from figures from the ADS article 
service dDemleitner et alj|200ll) . The selected lines and the result- 
ing AV are listed i n Tab les 1-2. 

ILedoux et alj J2006J) found the slope of the metallicity- Al 7 re- 
lation to be 1.46 and no evidence was found for an evolution of this 
slope with redshift. We shall therefore at first conservatively adopt 
this slope. This forms our conservative assumption, but in order to 
certify that our conclusions do not depend on this assumption we 
shall also repeat the entire analysis with less conservative assump- 
tions below. 

The sample is plotted in the left panel of Fig. 1 where we also 
plot the best fit linear relation for a slope of 1.46. The errors on 
[M/H] are typically in the range 0.05 — 0.20 and it is seen that 
the errors on the measurements are far too small to account for the 
scatter. As a first step we shall now determine the internal scatter of 
the distribution. To do this we assume that the measurement errors 
on the metallicity (<7 m et) as well as the natural internal scatter of the 
relation (cr nat ) both follow normal distributions and that the two are 
statistically independent. This allows us to define a total predicted 
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Table 1. Metallicities and velocity widths of low-ionization lines for z < 1 . 5 DLAs (metallicity measurements are from various papers quoted in the table 
footnote) 



Quasar 


--abs 


log TV (Hi) 


[X/H] 


X 


AV 
(kms^ 1 ) 


Selected 
transition lines 1 


Refs. 2 


B 0105—008 


1.3710 


21.70 + 0.15 


-1.40 + 0.16 


Zn 


31 


CrllA2056 


a 


J 0256+01 10 


0.725 


20.70 + 0.20 


-0.11 + 0.20 


Zn 


299 


MglA2852 


b 


Q 0302-223 


1.0094 


20.36 + 0.11 


-0.56 + 0.12 


Zn 


61 


SillA1808 


c 


Q 0449- 1645 


1.0072 


20.98 ± 0.06 


-0.96 + 0.08 


Zn 


204 


ABA 1862 


d 


Q 0454+039 


0.8597 


20.69 ± 0.06 


-1.01 + 0.12 


Zn 


111 


FeIlA2260 


c 


Q 0933+733 


1.479 


21.62 ± 0.10 


-1.58 + 0.10 


Zn 


30 


FeIlA2249 


e 


Q 0948+433 


1.233 


21.62 ± 0.05 


-1.14 + 0.06 


Zn 


65 


FenA2260 


e 


J 1009+0713 


0.1140 


20.68 ± 0.10 


-0.62 + 0.16 


S 


56 


TillA3384 


f 


J 1107+0048 


0.740 


21.00 + 0.04 


-0.54 + 0.16 


Zn 


212 


MglA2852 


b 


Q 1354+258 


1.4200 


21.54 + 0.06 


-1.61 + 0.16 


Zn 


25 


FellA2260 


g 


J 1431+3952 


0.6018 


21.2 ± 0.1 


-0.80 + 0.21 


Zn 


58 


MMIA2576 


a 


J 1623+0718 


1.3357 


21.35 ± 0.10 


-1.07 + 0.13 


Zn 


59 


FenA2260 


a 


J 2328+0022 


0.652 


20.32 ± 0.06 


-0.49 + 0.16 


Zn 


97 


MglA2852 


b 


B 2355-106 


1.1723 


21.0 ± 0.1 


-0.87 + 0.20 


Zn 


131 


MglA2852 


a 



1 Transition lines used to determine t he v elocity widths of low -ion ization line profiles . 

2 Refe renc es: rallEllison et all 120121) ; rbl lPeroux et"aTJto06h : rcl lPettini et all feOOOl) ; rdl lPeroux et"afl 120081) ; fel lRao et aljfe005h : rfl lMeiring et"afl 
<201ll) : rgl lPettini et afNl999h . 

Table 2. Metallicities and velocity widths of low-ionization lines for z ~ 4 DLAs (metallicity measurements are from lRafelski et all l|2012j) based on element 
X with an upward correction of 0.3 dex for Fe) 



Quasar 


•2-abs 


log AT (Hi) 


[X/H] 




AV 


Selected 










X 


(kms^ 1 ) 


transition lines 1 



J0040-0915 


4.7394 


20.30 ± 0. 


15 


-1.40 + 0.17 


Fe+0.3 


67 


FellA1608 


J 0747+4434 


4.0196 


20.95 ± 0. 


15 


-2.28 ± 0.20 


Fe+0.3 


69 


FeiiA1608 


J0817+1351 


4.2584 


21.30 + 0. 


15 


-1.15 ± 0.15 


S 


89 


SHA1250 


J 0825+3544 


3.2073 


20.30 ± 0. 


10 


-1.68 + 0.16 


Fe+0.3 


30 


FellA1608 


J 0825+3544 


3.6567 


21.25 + 0. 


10 


-1.83 + 0.13 


Si 


89/138 2 


SillA1808/SillA1304 


J 1051+3107 


4.1392 


20.70 ± 0. 


20 


-1.99 + 0.21 


S 


87 


SillA1304 


J 1051+3545 


4.3498 


20.45 ± 0. 


10 


-1.88 + 0.10 


Si 


29 


SillA1808 


J 1051+3545 


4.8206 


20.35 ± 0. 


10 


-2.28 + 0.10 


Si 


39 


SillA1526 


J 1100+1122 


4.3947 


21.74 + 0. 


10 


-1.68 + 0.18 


Fe+0.3 


137 


NillA1709 


J 1200+4015 


3.2200 


20.85 ± 0. 


10 


-0.64 + 0.10 


S 


145 


NillA1370 


J 1200+4618 


4.4765 


20.50 + 0. 


15 


-1.38 ± 0.16 


Fe+0.3 


88/113 2 


FeIlA1611/SillA1304 


J 1201+2117 


3.7975 


21.35 + 0. 


15 


-0.75 ± 0.15 


Si 


490 


SillA1808 


J 1201+2117 


4.1578 


20.60 ± 0. 


15 


-2.38 + 0.15 


Si 


25 


SillA1526 


J 1202+3235 


4.7955 


21.10 + 0. 


15 


-2.36 + 0.16 


Fe+0.3 


18 


Cll*A1335 


J 1202+3235 


5.0647 


20.30 + 0. 


15 


-2.66 ± 0.16 


Si 


30 


OIA1302 


J 1304+1202 


2.9131 


20.55 + 0. 


15 


-1.65 ± 0.16 


S 


64 


FeIlA1608 


J 1304+1202 


2.9289 


20.30 + 0. 


15 


-1.54 + 0.16 


S 


22 


FenA1608 


J 1353+5328 


2.8349 


20.80 ± 0. 


10 


-1.38 + 0.10 


S 


43 


SIIA1253 


J 1438+4314 


4.3990 


20.89 ± 0. 


15 


-1.31 ± 0.15 


s 


91 


FenA1608 


J 1541+3153 


2.4435 


20.95 ± 0. 


10 


-1.49 + 0.11 


Si 


47 


SinA1808 


J 1607+1604 


4.4741 


20.30 ± 0. 


15 


-1.71 ± 0.15 


Si 


37 


SillA1304 


J 1654+2227 


4.0022 


20.60 ± 0. 


15 


-1.66 + 0.16 


Fe+0.3 


22 


FenA1608 



1 Transition lines used to determine the velocity widths of low-ionization line profiles. 

2 For this system no line fulfilled the optimal criteria fully so we measured both a very weak line and a strong (slightly saturated) line. We list both lines and 
both measurements but we use the mean of the two. 



a 

O-tot(i) = (OmotM + Onat) 1/2 (1) 

for each observation, "i". For a given slope (o?o) and zero point (zp) 
of the metallicity-mass relation 

[M/H] = a log(AV) + zp (2) 



we can now assign a goodness of fit as the sum of the squares of 
the deviations of the fit to the data 

no 

C'L = ^2(([M/H](i) - a Iog(AV(t)) - zp)/a tot (i)) 2 )/dof 

(3) 

where the degrees of freedom (dof) in this case is 108. Cl 

of lS 

equivalent to the \ 2 P er degree of freedom of the distribution if 
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Table 3. DLA sample: average metallicities and velocity widths of low-ionisation line profiles * 



Quasar 




Other name 


2cm 


^•abs 


log/V(Hl) 


[X/H] 




AV 


Selected 


Refs. 2 
















X 


(tar 1 ) 


transition lines 1 




Q0151- 


r048 




1.93 


1.934 


20.34 ± 0.02 


-1.93 + 0.04 


Si 


51 


SillA1304 


a, b 


Q0918- 


-1636 


SDSS J 091826.16+163609.0 


3.07 


2.583 


20.96 ± 0.05 


-0.12 + 0.05 


Zn 


295 


SillA1808 


c 


Q0918- 


hi 636 


SDSS J 091826.16+163609.0 


3.07 


2.412 


21.26 + 0.06 


-0.55 + 0.16 


Zn 


352 


Si II A 1808 


d 


Q2222- 


-0946 




2.93 


2.345 


20.65 + 0.05 


-0.46 + 0.07 


Zn 


185 


SillA1808 


e 



Transition lines used to determine the velocity widths of lo w-i onization line profile s. 
2 References: faliZafar et al.N201 ll); fbl lEllison et aL1j2012l): MlFvnbo et aT] 1201 ll); rdllThorsenl fepi ll); fellFynbo et all 12010). 



(j nat is correctly chosen. We can therefore use eq(3) to select all 
pairs of (zp, <r nat ) for which C 2 of = 1, and we then choose the pair 
for which a nat has its minimum value as the globally best fit: (qo, 
zp, cr na t) = (1.46, -4.1 1 1, 0.4181). In the top right panel of Fig. 1 
we plot the residuals after subtraction of this fit. 



2.3 Redshift evolution 

By splitting their s ample into a high and a low redshift sample, 
iLedoux et al.l d2006t) already pointed out that while the slope of the 
correlation did not appear to change, there was a shift of the cor- 
relation. Here we seek a more general formulation of this redshift 
evolution. 

Inverting eq(2) we can interpret each absorber of the sample 
as an independent measurement of the zero point 



zp(i) = [M/H](i) - a log(AV(i)) 



(4) 



with an uncertainty of Otot(i). This allows us to check if we see 
an evolution of zp as a function of redshift. The function describes 
directly the evolution (shift) of the mass-metallicity relation at a 
given intercept (mass) and we have chosen to normalize the zero 
point to a mass of 100 km/s, which is a typical mass in our sample, 
and to name it [M/H] 100 k m / s (z). This function represents the 
mean metallicity of a galaxy with a mass of 100 km/s as a function 
of redshift. 

In Fig. 2 we show [M/H] 100 km/s( z ) plotted versus redshift, 
i.e. we plot the evolution of the mean metallicity of the cold gas in 
a galaxy of a mass of 100 km/s from redshift 5 to 0. The horizontal 
bars in Fig. 2 show the z range of each bin, while the vertical bars 
show the propagated uncertainty of the [M/H] 100 fc m / s (z). 

From Fig. 2 we see clear evidence for a steady and continuous 
incline of [Al/H] 100 km/s( z ) from redshift 2.6 to 0. An incline of 
M(z) means that the metallicity-mass relation changes such that 
for any given mass the galaxy becomes more metal rich at later 
times. At earlier times than redshift 2.6 there is no evidence for 
any evolution and it appears that the mass-metallicity relation was 
constant from z — 5 to z = 2.6. 

We will take this apparent "two phase" nature of the early uni- 
verse as our working hypothesis and proceed to fit the evolution 
with a constant at high redshifts and the late incline with a linear 
function. The transition between the two then happens at the tran- 
sition redshift z tr an and we get 



zp(z) 



Nearly (cOUSt) for Z > Ztran, 
CHate 2 + (Slate for Z < 2 t ran 



(5) 



where the constant value at high redshift is defined by the require- 
ment that zp(z) is continuous at zt ran , i.e. /3 e ariy is not a free pa- 
rameter. 
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Figure 2. Redshift evolution of the zero point of the DLA mass-metallicity 
relation. We have normalized the function to a gravitational well mass of 
100 km/s. The red line is a best fit to the data of the form given in eq(5). 
There is no evidence for evolution before redshift 2.6, and then a sharp 
break at z = 2.6 with an onset of strong evolution. The evolution is such 
that the metallicity of a galaxy of a given mass grows by 0.35 ± 0.07 dex 
per unit redshift. The fit was computed as a fit to all 110 data points, not to 
the binned data. 



Table 4. Fitted evolution parameters for ag 
defined in eq(5). 



1.46. The parameters are 



ttlate 


-0.35 ± 0.07 


Plate 


-3.33 + 0.16 


■2-tran 


2.62 + 0.20 


/^carly 


-4.25 ± 0.05 



We do not fit the binned data in Fig. 2, rather we proceed using 
the method described in Sect. 2.2 fitting simultaneously the param- 
eters in eq(2) and eq(5) for all 110 systems. In this way we find 
2tran =2.618 and cr na t=0.3759 for a =1.46 and 106 dof. In Table 4 
we provide all the parameters for the fit and in Fig. 2 this best fit is 
shown as red dashed lines. In Fig. 1 (lower panel to the right) we 
plot the residuals (res [m/h]) after subtraction of this fit. Comparing 
the upper and lower right frame it is easy to notice the decreased 
scatter of the residuals (from 0.4181 to 0.3759). 
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Figure 3. Same as Fig. 1 except that we here fit for the slope and find a best 
fit for a = 1.12. 

2.4 Testing the importance of the assumed an 

In the above analysis we started by fixing t he slope of the mass 
metallicity relation to a?o=1.46 as found by iLedoux et all {2006 ) 
who used a bisector fit. We agree that because of the large scat- 
ter the bisector is the best way of determining the slope, but we 
wish to test if another assumption would significantly change our 
result regarding the redshift evolution. We therefore now relax this 
constraint and repeat the analysis treating qq as a free parameter 
with an additional minimization fit. We find that the best fit in the 
form of eq(2) is obtained with a slope of Qo = 1.12 as shown in 
Fig. 3. The resulting fit to the redshift evolution is however almost 
unchanged. The transition redshift is now found at zt ran =2.622 (a 
change of only 0.003) and the redshift evolution slope is 0.37 dex 
per unit redshift (a change of 0.02). Both values are well within the 
la errors of the original fit. 

We conclude from this that the specific assumption of the 
slope of the underlying mass metallicity relation has almost no ef- 
fect on the results regarding the redshift evolution of the relation 
and that the result displayed in Fig. 2 is very robust. In hindsight we 
can now understand the riddle which for almost 2 decades was per- 
plexing cosmologists: "Why do we not see clear evid ence for cos- 
mic c hemical evolution in high redshift DLAs?" (e.g., Pett ini et"al] 
1 1994) . 

2.5 Comparison with the redshift evolution of the galaxy 
stellar mass-metallicity relation 

Maio lino et all (2008) provided galaxy mass-metallicity relations 
back to a redshift of 3.5 as a function of galaxy stellar mass (their 
eq. 2). We have computed a grid of those relations and overplot 
them on our DLA m easurements in Fig. 4. It can be seen that 
iMaiolino et alj {2008) produces relations of metallicity as a func- 
tion of redshift of increasing slopes with decreasing stellar masses. 
In the mass range log(M*/MQ) = 8.5 — 9.0, it is apparent that 
the slope matches the evolution of DLAs below z ~ 2.6 almost 




12 3 4 

Z abs 



Figure 4. Here we plot the reds h ift evo lution of the galaxy mass-metallicity 
relations from IMaiolino et all | |2008|) for different stellar masses (full 
curves) and the DLA data points and fit (dashed line) from Fig. 2. There 
is a remarkable agreement between the slope of the redshift evolution of the 
DLA galaxies and that of low-mass galaxies from Maiolino's study. 

perfectly. The fact that the evolution of the zero point of the DLA 
velocity-metallicity correlation with redshift matches the galaxy 
stellar mass-metallicity evolution so well strongly supports the idea 
that we are witnessing in the DLA data a mass-metallicity rel ation 
akin to galaxies as initially suggested bv lLedoux et alj feoOrj) fur- 
ther. 

It is still a matter of debate how exactly emission-line based 
metallicities relate to the metallicity of the absorbing gas. It is 
usually thought that the DLA gas can have lower metallicities 
than t hose determined from emission line fluxes (e.g. IPeroux et al.1 
120121) . However, there are uncertainties in the recip es used to mea- 
sure o xygen abundances from emission lines (e.g.. iKudritzki et all 
l2012h and in the choice of a suitable element to derive absorption- 
line metallicities. Other works have indeed found consistency be- 
tween the two measurements (e.g.. iBowen et alj|2005l : iKr ogager 
120130 . Here, we conservatively assume that the DLA gas either has 
the same metallicity or a metallicity up to 0.5 dex lower than that 
derived from emission lines. Any effects related to a-element over- 
abundance will also be covered by this assumption. Given the re- 
markable similarity of the slopes of the relations shown in Fig. [4] 
we then estimate that the "typical" DLA in our sample exhibiting 
AV = 100 km s -1 has a stellar mass log(Af*/Af0) of the order 
of 8.5. This is in agreement with the early result that due to the 
effect of gas cross-section selection even the most massive DLA 
gal axies only correspo nd to the least massive Lyman-Break galax- 
ies dFvnboetal.il 19991) . As a matter of fact, there is only a small 
range of overlapping metallic ity between our DLA sample and the 
galaxy sample in the study of Maioli no et alj j2008h . 

In order to generalize the comparison, we then determined the 
best fit to a range of models in the very lowest range of masses cov- 
ered by Maiolino's data, which constitutes the small overlap region 
with the DLA galaxies. In practice we created a finer grid of the 
Maiolino models than shown in Fig. 4, we then slid our data points 
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Figure 5. Offsets (A[M/H] 100km ^ s ) required to fit the data points in Fig. 4 
to a grid of stellar masses using the Maiolino models in the region of over- 
lap. The offsets are well fitted by a linear relation with slope 1.76 -1 . 

in Fig. 4 up to match each model defined by a given stellar mass 
and determined the offset (A[M/H] 10 ok m /s) required. The offsets 
are shown in Fig. 5 as a function of stellar mass and it is seen that 
in the mass range considered the offsets are well fitted by a linear 
relation with slope 1.76 - . 

The correction in metallicity when comparing absorption to 
emission is not yet well known so we allow for this as a free param- 
eter, C[ M /H]> which must be added to the absorption-line metallic- 
ity to make it consistent with the emission-line metallicity. Finally, 
in the redshift range considered we found an evolution with redshift 
such that A[M/H]= — 0.35Az for a given mass. We can now write 
up the relation: 

Iog(M»/M Q ) = 1.76([M/H] + G [M/H] + 0.352 + 5.04) (6) 

between DLA metallicity and stellar mass of the DLA host galaxy. 
Here 1.76 on the right side comes from the slope of the relation in 
Fig. 5, 0.352 from the redshift evolution (aiate in eq(5)) and 5.04 
is the empirically determined zero point of the relation. 

A number of important caveats should be noted before using 
this relation. 

• There is a large scatter in this relation (0.38 dex in [M / H]) so 
it is only useful as a statistical tool on large samples. Application 
to single systems will probably not be meaningful. 

• The relation is found from, and only known to be correct in, 
the small range of overlap between the two samples. 

• The relation is formally valid back to z — 2.6 and the fit 
is excellent between z = 2.6 and z — 0. At higher redshifts, 
flux-limited galaxy studies only sample the most massive galax- 
ies, while the DLA systems sample halos of roughly the same mass 
at all redshifts. There is therefore a huge gap between the two sam- 
ples at high redshifts and we have not tried to extract a similar rela- 
tion past the break at z — 2.62. The results provided in this paper 
suggest that the relation remains constant at higher redshifts (i.e., 
2 should be replaced by 2.62 in Eq. [6}, but this has not yet been 
verified by high-redshift galaxy studies. 



Pontz en et al. U20081) were able to confirm that the relation between 
AV and metallicity reported by iLedoux et al. (2006) is indeed a 
mass-metallicity relation. A relation between AV and virial mass 
of the DM halo is clearly seen in their simulations but the relation 
has a significant scatter. They also tend to underestimate the ve- 
locity widths of absorption-line profiles (their fig. 8). Taking this 
into account, a typical DLA with width 100 km s _1 would corre- 
spond to a virial mass log(M v i r /Mq) of the order of 10. This leads 
to a mass-to-light ratio of ~ 30 for a typical DLA galaxy at red- 
shifts below z ~ 2.6. We note that significantly larger virial masses 
have been inferred for DLA host halos from clustering arguments 
dFont-Ribera et alj|2012h . 



3 DISCUSSION AND CONCLUSIONS 

3.1 Transition redshift of the mass metallicity relation 

We confirm the reported mass metallicity relation for high redshift 
DLAs and w e confirm that there is indeed a redshift evolution of 
this r elation l lMeiller et ai]|2004 IL"edoux et al.ll2006t IPontzen et al.1 
2008) but only in the later phase of cosmic evolution. At high red- 
shifts the mass metallicity relation is constant while there is a break 
at 2 = 2.6±0.2 with a sudden onset of redshift evolution with a 
slope of 0.35 ± 0.07 dex per unit redshift. 

Previously similar sharp transitions i n the properties of Lyq 
emitters identified as ULIRGs and QSOs dNilsson & M0llerll2OQ9l : 
iBongiovanni et al.l I2OI0I ; iNilsson & Mailed 1201 lb an d in the de- 
tectio n rate of highly reddened quasar host galaxies (Fvnbo et al. 
120131) at redshifts around 2 = 2.5 — 2.6 have been reported. The 
close coincidence of the redshifts of those transitions, which also 
precisely coincides with th e peak of the star formation rate density 
( Hop kins & Beacoml F2006) makes it plausible that they must all be 
related in some way. Lya emitters only form a sub-set of the popu- 
lation of high redshift galax ies. DLAs however trace the gas both in 
Lya emitting galaxies (e.g. lM0ller& Warren] dl993h ) and in galax- 
ies wi th no Lya emission (e.g. Fvnbo et al] d201 lh : FBouche et"al] 
(2012)) and in fac t most DLA galaxies are non-detected in Lya 
(M0ller et al. 2004). It is therefore now evident that the transition 
seen for Lya emitters is not an artefact of the Lya emission selec- 
tion method. A recent detail ed study of the mass m etallicity relation 
of galaxies seen in emission Ma nnucci et al. I feoioh showed that the 
relation in fact forms a plane with the star formation rate as the third 
parameter. Their description forms a perfect fit in the redshift range 
2 = — 2.5, but breaks completely down at z=3 as seen in their 
Figure 4, right panel. 



3.2 Slope of the redshift evolution 

The redshift evolution below 2 = 2.6 is such that for a given virial 
halo mass the gas has higher metallicity at lower redshifts than at 
higher red shifts. This is in agr eem ent with results for g alaxies re- 
ported by Savagl io et al. U2005I) and Maiolino et al. (2008). We find 
the slope of the evolution to be 0.35 ± 0.07 dex per unit redshift 
in agreement with the slope found in the low stellar mass range for 
similar studies of galaxies. 



2.6 Virial masses 

We now ask the question what the virial mass (Af v i r ) of the 
dark matter halo of a typical DLA galaxy is. Via identification 
of DLAs in dark matter halos of their cosmological simulations, 



3.3 The stellar masses of DLAs 

Identifying the galaxies that harbour the DLAs at high redshifts 
has been a holy grail for DLA studies ever since the first DLA was 
discovered. Attempts to do this on a case by case basis have been 
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mostly unsuccessful (e.g. JSmith et alj 19891) and still only a coupl e 
of handful detections have been reported dKrogageretaLll2012h . 
In this paper we were able to tie a direct link between a flux se- 
lected galaxy sample and our DLA sample using the evolution of 
the mass-metallicity relation to link them, and to provide an expres- 
sion for the stellar mass of DLAs as a function of their metallicity 
and redshift. 



3.4 A plausible cause of the transition: Infall starvation 

The conclusions from a study of the redshift evolution of Lya emit- 
ters dNilsson et al]|2009h were that at redshifts above ~ 3 the ob- 
jects were all star forming, with little or no dust, and with young 
blue populations, while at redshifts lower than ~ 2.5 a signifi- 
cant fraction of red, dusty, mixed population objects, even ULIRGs, 
were found. The simplest way to explain those observations as well 
as the new results reported here is a scenario in which primor- 
dial, or near-primor dial, gas is constantly accreted onto the galax- 
ies at high redshif ts jB arkana & L oeb 2003; IWeidinger et al. 2005; 
iDekel et alJ2 009). If the amount of gas available for infall is larger 
than a galaxy can accept and process then they will all accrete at 
maximum capacity. As the gas is converted into stars and enriches 
the galaxy with metals, this sets up an unchanging universal rela- 
tion between the mass and the metallicity of a galaxy. With respect 
to this relation one might refer to this early phase as a steady state 
of galaxy evolution. 

In order to change the universal mass-metallicity relation, 
some other parameter has to change first and the only available 
parameter is the infall rate. We therefore propose that the trigger 
to start the evolution of the mass-metallicity relation is a drop in 
the infall rate of primordial gas at Ztran = 2.6 below the point of 
balance for the steady state to be maintained. 

This can be accomplished in two ways. Either the supply 
of metal poor gas simply drops far below that of metal en- 
riched gas, or at the peak of the star formation rate at z ~ 2.6 
(Hopkins & Beacon] [2006) powerful galactic winds create highly 
ionized halos which render it difficult for the infall to continue at 
the same high le vel. This effect is k nown as accretion quenching 
and was used bv lBouche et al. I d2010h who found that it was a nec- 
essary ingredient in their model in order to reproduce the observed 
star formation rate-mass and Tully-Fisher relations. The primary 
mode of galaxy stellar mass growth then shifts from ga s accretion 
to mer ging. There is some theoretical evidence for this. lOser et all 
reported that their simulations showed that there is a transi- 
tion between redshifts of 3 and 2 where the galaxy growth changes 
from "in situ" formation of stars from accreted cold gas towards 
accretion of stars which were formed in satellites earlier. The new 
results we have presented here are observational support for this 
model as wel l as fo r the accretion efficiency controlled model by 
iBouche et ailfeOld) . There are similar obvious links to the sudden 
change in Lya emission line selected samples at jz tra n = 2.52 and 
to the setup of the tight 3-way relation between mas s, star formation 
and m etallicity (FMR) between z — 3.3 and 2.5. iMannucci et al.1 
l2oToh . In the latter work it was found that infall rates must have 
been ~ 10 3 times higher at * = 3.3 than locally. 



3.5 Future work and outlook 

The dominating source of errors is the internal scatter, not the mea- 
surement errors. Using a larger sample or a sample with smaller 
errors will therefore not improve the results much. The best way to 



obtain a significant improvement is to identify and understand the 
source of the scatter. 

N-body simulations of the galaxy growth via "in situ" forma- 
tion of stars versus accretion of satellites including gas metallicity 
predictions would be an important tool to test if the proposed model 
is correct. 

In order to be able to test if the break at z = 2.6 is also 
present in galaxy emission selected samples those would need to 
reach much lower stellar mass at redshifts between 4 and 5 to make 
it possible to tie them together with DLA samples. This will likely 
not be possible until E-ELT class telescopes become available. 
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